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Abstract. Unpulsed, high-energy emission from pulsar binaries can be attributed to the interaction 
of a pulsar wind with that of a companion star. At the shock between the outflows, particles 
carried away from the pulsar magnetosphere are accelerated and radiate both in synchrotron and 
inverse Compton processes. This emission constitutes a significant fraction of the pulsar spin-down 
luminosity. It is not clear however, how the highly magnetized pulsar wind could convert its mainly 
electromagnetic energy into the particles with such high efficiency. Here we investigate a scenario 
in which a pulsar striped wind converts into a strong electromagnetic wave before reaching the 
shock. This mode can be thought of as a shock precursor that is able to accelerate particles to 
ultrarelativistic energies at the expense of the electromagnetic energy it carries. Radiation of the 
particles leads to damping of the wave. The efficiency of this process depends on the physical 
conditions imposed by the external medium. Two regimes can be distinguished: a high density one, 
where the EM wave cannot be launched at all and the current sheets in the striped wind are first 
compressed by an MHD shock and subsequently dissipate by reconnection, and a low density one, 
where the wind can first convert into an electromagnetic wave in the shock precursor, which then 
damps and merges into the surroundings. Shocks in binary systems can transit from one regime to 
another according to binary phase. We discuss possible observational implications for these objects. 
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INTRODUCTION: THE cr-PROBLEM 

A wind, outflowing from a pulsar magnetosphere, is a mixture of plasma, produced 
close to the star, and the strong electromagnetic fields that are anchored in its surface 
and twisted beyond the light cylinder = cj CO (CO is the angular velocity of the star). 
This wind carries the entire pulsar spin down power L, expressed via the dimensionless 
parameter «l = (e 2 L/m 2 c 5 ) 1//2 , an d me mass loading of the wind is given by \l = 
Ol/Ak, where K is the multiplicity of pulsar cascades. The energetics of the outflow 
is completely dominated by the Poynting flux, so that the magnetization parameter, 
the ratio of the energy flux carried by the fields to that carried by the particles, is 
large o ^> 1. Within the ideal MHD description of the wind, a large initial value of 
o remains so along the entire flow, up to the point where the wind impacts on the 
surroundings and terminates in the standing shock. However, the observed emission, 
attributed to the synchrotron processes downstream of this shock, suggests that most of 
the energy is deposited in the particles. The mechanism that transfers the energy stored 
in the fields to the plasma is not clear, and called the "a-problem". A solution must 
be envisaged beyond the domain of validity of the ideal MHD approach. In particular, 
as the wind propagates radially outwards, the plasma density decreases, and beyond 
the critical radius r c = a^r^/ /A, it is low enough to enable the propagation of strong 
electromagnetic (EM) waves. Such an underdense medium, in contrast to the MHD 



plasma, cannot screen out the wave displacement currents, and when they take over, 
mode conversion between the MHD wind and the EM mode can occur. This transition 
can be regarded as a consequence of the boundary conditions, and the new EM wave 
mimics a precursor to the shock, launched to adjust the flow to the surroundings. This 
mode is able to efficiently transfer the EM energy to the plasma, providing a solution to 
the a -problem. 

ELECTROMAGNETIC SHOCK PRECURSORS 

Properties and radial propagation of large-amplitude waves. Propagation of strong 
EM waves is described within the two-fluid model of a plasma that is coupled to the 
EM fields via equations of motion and Maxwell equations. Under pulsar conditions the 
plasma is positronic, which implies that the electric vector of EM waves is purely trans- 
verse. The e ± fluids move with equal momentum pn = yvu/c in the propagation direc- 
tion, but have oppositely directed, equal amplitude oscillations in transverse momentum 
p±. In the lab. frame these waves have superluminal phase velocity /3 > 1, but sublu- 
minal group speed /3* = 1/j8 < 1 that, in general, does not coincide with the parallel 
component of the fluid 3-velocity vy/c (i.e., the particles stream through the wave). 

Using perturbation analysis in the small parameter £ = r^/r <C 1 (short wavelength 
approximation; r is the distance at which the wave is probed), we expand all the quan- 
tities in order to obtain a set of equations that, to the lowest order, describe the plane 
wave solution, dependent only on the wave phase. To first order they govern the slow 
radial evolution of the phase-averaged, plane- wave quantities. The radial propagation is 
given by (1) the continuity equation, (2) energy conservation, and (3) the evolution of 
the radial momentum flux. In contrast to the MHD wave, the radial momentum flux is 
not conserved in spherically expanding EM modes, because the flow exerts a pressure in 
the lateral direction, and the energy division between the parallel and these perpendicu- 
lar degrees of freedom changes. In addition to the conserved particle and energy fluxes, 
it can be shown [1] that the phase- averaged Lorentz factor of the particles, measured 
in the laboratory frame, ("Hah)* is an integral of motion for both circularly and linearly 
polarized modes. 

The EM mode at the conversion radius is determined from the jump conditions be- 
tween the MHD and the EM wave, which carry the same particle, energy and radial mo- 
mentum fluxes [2, 3]. Fig. 1 shows the Lorentz factor of a strong EM wave y*, obtained 
from the jump conditions (dashed curves), and its radial evolution (solid curves) for dif- 
ferent launching radii. There are two solutions of the jump conditions that describe two 
possible EM modes: a free-escape mode (higher branch) and a confined mode (lower 
branch). Only the latter decelerates at large distances, and, since the wind is terminated 
in the (roughly) standing shock, in the following we concentrate on this mode. 

The unique determination of a wave at launch requires the knowledge of four quanti- 
ties: the conversion radius Rq, initial group speed /3*o, and initial particle momenta pp, 
Plq- The jump conditions define three of them, leading to the dotted curves in Fig. 1. 
However, the evolution equations of the EM mode imply that its ram pressure tends to 
a constant value at large radius, which, in fact, should be equal to the external pressure 
Pext- Since only the wave launched at the correct radius has the correct asymptotic pres- 
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FIGURE 1. Dashed: The Lorentz factor of a circularly polarized EM wave, corresponding to an MHD 
wind with fi = 10 4 , a = 100, a L = 3.4 x 10 10 ; against conversion radius. Solid: the radial evolution of an 
EM wave launched at a point on the dashed curve (marked by a dot). Radius is normalized R = r^i(o/(caT_). 

sure, the fourth quantity, Rq, is determined by the matching procedure. Given that ('tfab) 
is conserved, one can connect the asymptotic wave pressure with the initial wave pa- 
rameters. Therefore, a unique, stationary precursor-solution is determined by the MHD 
wave parameters o, fx, a^, and also p ext , unless p ext is larger than the ram pressure at 
r = r c , p(r c ) = p c , in which case no such precursor can be formed. 

Radiation damping of electromagnetic waves. Self-consistent EM waves are prone 
to damping by nonlinear inverse Compton (NIC) scattering (i.e., radiation reaction of 
the particles accelerated in the wave fields) [4], as well as inverse Compton scattering 
of ambient low energy photons. Even though the flow energy is dissipated to the photon 
field, the particles themselves gain more energy from the fields than they lose in radi- 
ation, and both the wave and the particle streaming slow down due to the interaction. 
Thus, this effect resembles that caused by radial expansion. When, by either of these 
mechanisms, the streaming of the particles through the wave vanishes, strong waves be- 
come very unstable to small density perturbations in the direction of motion [5, 6, 7]. 
Parametric instabilities then set in and destroy the wave in only one wave period. This 
point of vanishing streaming is the location of the termination shock, beyond which the 
flow energy is thermalized. 

IMPLICATIONS FOR THE EMISSION IN BINARIES 

For specified values of the magneto spheric and wind parameters, the shock structure is 
determined only by the pressure conditions in the external medium. We identify two 
regimes: 

one with high external pressure p ext > p c - superluminal mode cannot be launched, 
because the conversion radius, required by the external conditions, falls below the 
critical radius; the current sheets in the pulsar wind are first compressed by an MHD 
shock and subsequently dissipate by reconnection [8], 



one with lower external pressure p ex t < p c - the EM precursor can be launched; in 
this case two possibilities arise: 

(1) Pext < Pc - radiation damping of a wave is efficient and determines the shock 
location at the point of vanishing streaming, 

(2) Pext <^ Pc - radiation damping ineffective; the particle streaming decreases due 
to, e.g., radial expansion, until at large distances it drops sufficiently, the intrinsic 
instabilities set in and disrupt the wave, forming a shock. 

In binaries, the high-pressure environment is provided by a wind of a companion 
star. If in the periastron p ex t < p c , so that the precursor can be launched, it has a large 
amplitude, and therefore the NIC process very quickly decreases the wave streaming. In 
this case the shock is formed without significant emission from the precursor. Most 
of the flow energy can be deposited in the shock- accelerated particles and one can 
expect enhanced emission from the shock itself. Another possibility is that at periastron 
Pext > Pc, and the precursor can be launched only in a more distant part of the orbit. At 
larger orbital separation, the damping is less abrupt and the emission from the precursor 
may constitute a significant fraction of the radiative output of the system. 

In the binary PSR B 1259-63 the pulsar is on a very elongated orbit around a compan- 
ion star and the external conditions for the pulsar wind change with the orbital phase. 
Pext can be determined from the standard model of Be-star winds, e.g., [9]. The precur- 
sor can exist if the value of /i is high enough to support the EM wave solutions required 
by Pext- For this system, if ji ~ 10 5 the precursor can be launched in every point of the 
orbit. NIC emission in the periastron would be in the X-ray band, and in the optical 
band in the apastron - possibly hidden beneath the luminous stellar emission. However, 
if /i ~ 10 4 , the shock regime may be switched, with the MHD structure close to the pe- 
riastron, and an EM, non-MHD precursor further out. If the regime switching occurs 30 
days after the periastron passage, when the binary members are 3 x 10 13 cm apart, the 
critical radius equivalent to this distance would imply jU = 2 x 10 4 . This is equivalent 
to the multiplicity K ~ 10 5 , about an order of magnitude smaller than the one inferred 
for the Crab pulsar. The emission from the precursor would be in the optical band, and 
the IC emission associated with the particles upscattering the stellar photons, in a few 
hundred MeV - consistent with the high energy flare observed by Fermi-LAT [10]. 
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